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We report results of autocorrelation measurements in pure 5(7(3) lattice gauge theory. The computations are 
performed on the CONVEX SPP1200 parallel platform within the CANOPY programming environment. The focus 
of our analysis is on typical autocorrelation times and optimization of the mixing ratio between overrelaxation 
and pseudo-heatbath sweeps for generating gauge field configurations. We study second order tadpole-improved 
approximations of the Wilson action in the gluon sector, which offers the advantage of working on smaller lattices 
(8 3 x 16 and 6 3 x 12 - 30). We also make use of anisotropic lattices, with temporal lattice spacing smaller 
than the spatial spacing, which prove useful for calculating noisy correlation functions with large spatial lattice 
discretization (of the order of 0.4 fm). 



NOTATION 

• s - spatial 

• t - temporal 

• Tint ~ integrated autocorrelation time 

• £ - correlation length 

• a s - spatial lattice spacing 

• at - temporal lattice spacing 

• p = at/a s - anisotropy ratio 

• u s - tadpole-improved mean spatial link 

• ut - tadpole-improved mean temporal link 

• S n d - lattice action without doublers 

• Pss - spatial plaquette 

• Pst - temporal plaquette 

• Rss - spatial 2x1 rectangle 

• Rstl - temporal 2x1 rectangle of type-I 

• Rst2 - temporal 2x1 rectangle of type-II 
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• Act - average gluon action 

• H-B - heatbath 

• Over. - overrelaxation 

• MC - Monte Carlo 

2. MOTIVATION 

Recently it has become apparent that coarse 
anisotropic lattices with p < 1 can be very use- 
ful in performing accurate Monte Carlo simula- 
tions of QCD at low computational cost. This 
is especially true when modeling heavy quantum 
states of QCD - like glueballs, for example. Be- 
cause of the exponential fall off of the signal, small 
at gives better resolution of the correlators at an 
early time step. On the other hand, a s should be 
kept relatively large because of critical slowing 
down. 

Since the anisotropy of the lattice breaks the 
Euclidean invariance of the continuum theory, it 
induces temporal £j at and spatial £' at correlation 
lengths which scale as 
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so that at <C a. 
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On the other hand, the autocorrelations in 
Monte Carlo updates are proportional to a power 
of the correlation length 

(2) 



oc 



(Aat 

V^trclevcnt op. 



2 



where theoretically n = 2 for local stochastic up- 
dates and n — for cluster/overrelaxation up- 
dates. 

In practice, different lattice operators will have 
very different autocorrelation times, and we ex- 
pect operators that couple strongly to £j at to have 
larger r int . Unfortunately "interesting operators" 
in lattice QCD, like those for glueballs, live in the 
spatial domain and scale with £j at which is large. 

This work tries to address the issue of the 
scaling behavior of different gluon operators on 
anisotropic lattices, and their relevance to the 
problem of MC algorithm optimization. 
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3. THE GLUON ACTION 

The doubler-free gluon action used in this 
study is given by 
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This is a Symanzik- improved lattice action jjj], 
which is accurate up to errors 0(a A sl a%) (classi- 
cally) . 

The tadpole-improvement scheme is 
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because of the assumption a t <C a s (||| . 

We examine anisotropic lattices with spatial 
lattice spacings a s from 0.25 to 0.4 fm and 
anisotropy p ranging from 1/2 to 1/5. The mean 
spatial tadpole-improvement factor u s varies from 
0.772 to 0.807. For every lattice considered 
we generate a corresponding Markov chain of 
4000 configurations (i.e. MC sweeps), which are 
viewed as Monte Carlo time series of appropriate 
lattice operators. 

On these time series we perform autocorrela- 
tion measurements on six lattice observables: Pss, 
Pst, Rss, Rstl, Rst2 and Act, which are shown in 
Fig. | 



Figure 1. Types of planar loops used in this study. 



4. AUTOCORRELATION FUNCTIONS 

In Fig. H and Fig. || we present comparison 
graphs of the autocorrelation functions for the 
simple plaquettes and the 2x1 rectangles. At 
small anisotropy of the lattice the autocorrelation 
times of the spatial and temporal gluon operators 
are small and almost indistinguishable. 

With increasing (inverse) anisotropy, p , the 
degeneracy lifts. The autocorrelation times in- 
crease for Pss and Rss, and remain small and 
degenerate for Pst, Rstl and Rst2 observables. 
This is to be expected, since Pss and Rss couple 
strongly to £| at . 

As can be observed from the graphs our auto- 
correlation functions distinguish well between the 
spatial and temporal modes on the lattice, and 
can be a good starting point for the optimization 
procedure. 

5. AUTOCORRELATION TIMES 

The integrated autocorrelation times for an ob- 
servable X are obtained using the method pro- 
posed by Sokal pLp}, namely 



Tint = 7: 
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Figure 2. Effects of anisotropy on the autocorre- 
lation functions for the plaquettes. 



Figure 3. Effects of anisotropy on the autocorre- 
lation functions for the 2x1 rectangles. 



with 



R(t) = 
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and m chosen so that T mt <C m <C n. The small- 
est value of m for which m/T m t > 4 has been 
chosen in a self-consistent manner. From here we 
derive an estimate for the error of Tint given by 
the formula 
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We call this criterion - "criterion m/4". 

We also double check our results against the 
procedure employed by the QCD-TARO Collab- 



oration B, wherein the autocorrelation time is 
defined as: 
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Tint = P (0) + 2j2p(t) 
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where N is determined so that r; n t is maximized, 
but N < 10% of the total sample and < 3Ti nt . 
Likewise, we call this criterion - "criterion r max " . 

We use hybrid overrelaxation updating scheme, 
which is considered to be the state-of-the-art al- 
gorithm for bosonic systems when no cluster al- 
gorithm is available 0. It simply consists in 
the mixing of Cabibbo-Marinari pseudo-heatbath 
and the Brown- Woch microcanonical overrelax- 
ation sweeps with a ratio 1 : n a MM- 
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6. ANISOTROPY-DEPENDENCE OFr int 



7. CONCLUSIONS 



In Fig. |] we show our final results for the de- 
pendence of the integrated autocorrelation time 
on the anisotropy of the lattice for the six lattice 
observables. We find small variation of the am- 
plitude of the integrated autocorrelation time T in t 
for the temporal loops. The error bars for these 
observables are small. It is apparent that the 
modes are strongly suppressed because they effec- 
tively propagate in the spatial domain with char- 
acteristic correlation length £j, at , which is small. 
In contrast, the integrated autocorrelation times 
for the spatial loops increase as p _1 increases. 



m/4 - criterion 




1/3 1/4 
Anisotropy 
tau(max) - criterion 



1/6 




1/3 1/4 
Anisotropy 



Figure 4. Scaling behavior of r m t on anisotropy. 



This study is an initial attempt to perform ac- 
curate studies of QCD gluon update dynamics 
on lattices with large spatial lattice spacing and 
small temporal lattice spacing. 

This work is at an early stage of development. 
We have found clear evidence that the spatial 
and temporal link variables have different au- 
tocorrelation behavior, and that this difference 
becomes more extreme as the anisotropy is in- 
creased. We hope to look for an optimal up- 
date scheme for anisotropic simulations. This 
will be accomplished by exploiting the freedom 
we have on such lattices to vary pseudo-hcatbath 
and overrelaxation updates on the spatial and 
temporal links independently. Since the action 
has explicitly broken Euclidean symmetry, we can 
expect that an optimal update scheme will in- 
volve different update methods for spatial and 
temporal links. We anticipate the optimal update 
method will involve overrelaxation of the spatial 
links (which have been shown here to have the 
longest autocorrelation times) with heatbath up- 
dates of temporal links to ensure the overall up- 
date is ergodic. 
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